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We report a measurement of the amplitude ratio rs of B° D°K*° and B° D°K*° decays 
with a Dalitz analysis of D Kg7r + 7r~ decays, for the first time using a model-independent method. 
We set an upper limit rs < 0.87 at the 68 % confidence level, using the full data sample of 711 fb -1 
corresponding to 772 x 10 6 BB pairs collected at the Y(4 S) resonance with the Belle detector at 
the KEKB e + e _ collider. This result is obtained from observables X- — +0.4to.6io’.i =*= 0-0? V- — 
-0.6l°;®t°;J=b0.1,x+ = +0A±° o ]lt° o °±0A and +0.3tg;ftg;°±0.1, where x± = r s cos(fe±</> 3 ), 
y± = rssm(5s =b (^ 3 ) and <f 3 (5s) is the weak (strong) phase difference between B° D°K*° and 
B° -> D°K*°. 


PACS numbers: 
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INTRODUCTION 


Determination of parameters of the standard model (SM) plays an important role in the search for new physics. In 
the SM, the Cabibbo-Kobayashi-Maskawa (CKM) matrix (T] gives a successful description of current all measurements 
of CP violation. The CP-violating parameters </q, <p 2 and </>3 are the three angles of the most equilateral of the CKM 
unitarity triangles, of which </>3 = arg (— V u dV u b*/V c dV c b*) is the least accurately determined. In the usual quark- 
phase convention, where the complex phase is negligible in the CKM matrix elements other than V u b and Vtd [2], the 
measurement of (j )3 is equivalent to the extraction of the phase of V u b. To date, 0 3 measurements have been performed 
mainly with B meson decays into D^*)p(*) final states pMT2]. all of which exploit the interference between the P)(*)° 
and P)(*)° decaying into a common final state. In particular, Dalitz analyses of B ± -+ D^P^*^, D -+ K® s 7t~ ± tt~ 
provide the most precise determination of 0 3 . The Dalitz analysis technique for the measurement of </>3 was proposed 
in Ref. [13]. Belle reported the first 0 3 measurement with the model-independent Dalitz analysis technique in Ref. [13], 
which exploits a set of measured strong phases instead of relying on a D decay model into a three-body final state. 

In this paper, we present the first measurement of the amplitude ratio of B° -+ D°P*° and P° -+ D°P*° decays 
with a model-independent Dalitz analysis. We reconstruct B° -+ DP* 0 , with P*° -+ P + 7 r _ (Throughout the paper, 
charge-conjugate processes are implied; P*° refers to P*(892)° and D refers to either D° or D° when the D° flavor 
is untagged). Here, the flavor of the B meson is identified by the kaon charge. Neutral D mesons are reconstructed 
in the P^ 7 r + 7 r _ decay mode. The reconstructed final states are accessible through b -+ c and b -+ u processes via the 
diagrams shown in Fig. 
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In this analysis, we use the variables rs, k, and 5s to parameterize the strong dynamics of the decay. These 
parameters are defined as m 


r s - 


ke i5s 


r(B° -> D°K+n-) _ fdpA^Jp) 
T (B° -> D°K+ir~) f dpA^ c (p) ’ 

]/.f dpALJp) f CPCa (P) 


( 1 ) 

( 2 ) 


where the integration is over the B° —>• DK + ii~ Dalitz distribution region corresponding to the P*° resonance. Here, 
Ab^ c {Ab^u)(p) is the magnitude of the amplitude for the b —>> c (u) transition and S(p) is the relative strong phase, 
where the variable p indicates the position within the DK + tt~ Dalitz distribution. If the B° decay can be considered 
as a DP * 0 two-body decay, rs becomes the ratio of the amplitudes for b u and b —)► c and k becomes 1 . According 
to a simulation study using a Dalitz model based on the measurements in Ref. m, the value of k is 0.95 ±0.03 within 
the phase space of the DP * 0 resonance. The value of rs is expected to be around 0.4, which corresponds naively to 
I VubV* s I / I VcbV* s | but also depends on strong interaction effects. For rs, the best experimental value is reported 


by LHCb [171 as rs = 0.2407 


+0.055 


0.048 (different from zero by 2.7 a) from B° 


DP* 0 , D -+ P+P- 


K ± tt t decay. 


THE MODEL-INDEPENDENT DALITZ ANALYSIS TECHNIQUE 

The amplitude of the B° -+ DP* 0 , D -+ P^ 7 r + 7 r - decay is a superposition of the B° -+ D°P*° and B° -+ D°P*° 
amplitudes 


AB(rn\,m 2 _) = A ± rse^ Ss+ ^A, 


( 3 ) 
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where m\ and m 2 _ are the squared invariant masses of the K^tt^ and Kg7r~ combinations, respectively, A = 
A(m\,m 2 _) is the amplitude of the B° D°K *°, D° decay and A = A(m\,m 2 _) is the ampli¬ 
tude of the B° —>• D°K *°, D° —)> decay. In the case of CP conservation in the D decay, we have 

A(m+, m 2 _) = A(m 2 _, m 2 ^) as a CP transformation changes 7 ^ —>• 7 t t , thus m± m 2 1 . The Dalitz distribution 
density of the D decay from B° DK*° is given by 

Pb =| A b | 2 =| A + rse^s+^A | 2 = P + r 2 s P + 2 kVpP(x+C + y+S ), ( 4 ) 

where P = P(ra^_, m 2 _) =| A | 2 , P = P(ra^_, m 2 _) =| A | 2 , and 

x + = r s cos (S s + <fe), y+ = r s sin (S s + <fo). (5) 

The functions m?_) and m?_) are the cosine and sine of the strong-phase difference m 2 _) = 

arg A — arg A between the D° —» and D° amplitudes. Here, we have used the definition of 

k given in Eq. Q. The equations for the charge-conjugate mode P° —>> M*° are obtained with the substitution 
—0 3 —* 0 3 and dol; the corresponding parameters that depend on the P° decay amplitude are 


x- = r s cos (S s - 0 3 ), y- = r s sin (S s - 0 3 ). ( 6 ) 

If P, P, C, S and k are known, one can obtain (x+, y+) from B° and (x_, y_) from P° decays. Combining both P° 
and P° measurements, r*s, </> 3 and £ 5 - can be extracted. 

In the model-dependent analysis, one deals directly with the Dalitz distribution density and the functions C and 
S are obtained from a model based upon a fit to the D° — K^tt^tt - amplitude. On the other hand, in the model- 
independent approach [15] , where the assumption of a model for D° —>> K^tt^tt - decay is not necessary, the Dalitz 
plot is divided into 2A f bins symmetric under the exchange m 2 _ 0 m 2 ^. The bin index i ranges from —J\f to J\f 
(excluding 0); the exchange m 2 _ ^ m 2 + corresponds to the exchange i —i. The expected number of signal events 
in bin i of the Dalitz distribution of the D mesons from P° is 


N?=h B 


K±i + r 2 s K Ti + 2 ky/KiK-i(x±Ci ± y±st) , 


(7) 


where N + (A stands for the number of B°(B°) meson decays, 1 is the normalization constant and K +l is the 
number of events in the i th bin of a flavor-tagged D° —>• Kstt^tt - decays measured with a sample of inclusively 
reconstructed P* + — >► P° 7 T + decays. Equation ([7]) can be obtained by integrating Eq. Q over the i th bin region. 
Here, Ki oc J v |H| 2 dP, and V represents the Dalitz plane and Vi is the bin over which the integration is performed. 
The values of Ki are measured from a sample of flavor-tagged D° mesons obtained by reconstructing D * ± —)> P 7 r ± 
decays. The terms c$ and Si are the amplitude-weighted averages of the functions C and S over the bin: 


! Vi \A\\A\ CAD 


(8) 


The terms si are defined similarly with C substituted by S. The absence of CP violation in the D decay implies 
Ci = C-i and Si = —S-i. The values of q and Si can be measured using quantum-correlated D pairs produced at 
charm-factory experiments operating at the threshold of DD pair production. The CLEO Collaboration has reported 
Ci and Si values from CP tagged and flavor tagged DD events data, and this analysis is performed with the optimal 
binning in Refs. [El EDI, as shown in Fig. [ 2 } Given that q and Si are measured and Ki and k are known, Eq. 0 has 
only three free parameters (x, y, and hs) for each of B° and P°, and can be solved. We use the values of (c*, Si) for 
the “optimal D° —>> K^tt^tt - binning” reported in Table XVI of Ref. [20] . “the optimal binning” Ki values reported 
in Table II of Ref. [14], and k = 0.95 db 0.03 [IS] , We have neglected charm-mixing effects in D decays from both the 
B° —>• DK*° process and in the quantum-correlated DD production EH- 


EVENT RECONSTRUCTION AND SELECTION 

This analysis is based on a data sample that contains 711 fb -1 corresponding to 772 x 10 6 BB pairs, collected with 
the Belle detector at the KEKB asymmetric-energy e + e - (3.5 on 8 GeV) collider [ 22 ] operating at the Y(45) resonance. 
The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector, a 50-layer 
central drift chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement of 
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Phase bins 



FIG. 2: Binning of the optimal binning in Ref. 19, 2(1. 


time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return located outside of the 
coil is instrumented to detect K £ mesons and to identify muons. The detector is described in detail elsewhere [23] . 

We reconstruct B° — DK*° events with iT*° —> K + tt~ and D —>• Kg7r+7r~. The event selection described below 
is developed from studies of continuum data taken at center-of-mass energies just below the Y(4 S) resonance and 
Monte Carlo (MC) simulated events. 

The Kg candidates are identified using the output of a neural network. Inputs to the network for a pair of oppositely- 
charged pions are the invariant mass, 20 kinematic parameters and particle identification (PID) information from the 
ACC, TOF and the ionization energy loss in the CDC. The Kg selection has a simulated purity of 92.2% and an 
efficiency of 75.1%. Charged kaon and pion candidates are identified using PID information. The efficiency is 80- 
90% and the probability of misidentification is 6-10%, depending upon the momentum of hadrons and obtained 
using dedicated data control samples. We reconstruct neutral D mesons by combining a Kg candidate with a pair 
of oppositely-charged pion candidates. We require that the invariant mass be within ±15 MeV/c 2 (±3cr) of the 
nominal D° mass. K*° candidates are reconstructed from K + ii~ pairs. We require that the invariant mass be within 
±50 MeV/c 2 of the nominal iC*° mass. We combine D and iC*° candidates to form B° mesons. Candidate events 
are identified by the energy difference A E = Ei — E\> and the beam-constrained mass Mb c c 2 = \ c^2 i pi | 2 , 

where is the beam energy and pi and Ei are the momenta and energies, respectively, of the B° meson decay 
products in the e + e _ center-of-mass (CM) frame. We select events with 5.21 GeV/c 2 < M^ c < 5.29 GeV/c 2 and 
-0.10 GeV < A E < 0.15 GeV. 

Among other B decays, the most serious background is from B° decaying to the same final state as B° —)> DK*°. 
To suppress this background, we exclude candidates for which the invariant mass of the iC*°7r + system is within 
±4 MeV/c 2 of the nominal D + mass. This criterion leads to a negligible contamination from this mode and a relative 
loss of 0.6% in the signal efficiency. 

The large combinatorial background of true D° and random K + and tt~ combinations from the e + e - —>• cc process 
and other BB decays is reduced if D° candidates that are a decay product of D* + —)> D°7r + are eliminated. We use 
the mass difference AM between the [iC^7r + 7r - ]£)7r + and [Kg7r+7r~]]j systems for this purpose: if AM > 0.15 GeV/c 2 
for any additional 7r + candidate not used in the B candidate reconstruction, the event is retained. This requirement 
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1 

Fisher discriminants based on modified Fox-Wolfram moments [24]. 

2 

The angle in the CM frame between the thrust axes of the B decay and that of remaining particles. 

3 

The signed difference of the vertices between the B candidate and the remaining charged tracks. 

4 

The distance of closest approach between the trajectories of the K* and D candidates. 

5 

The expected flavor dilution factor described in Ref. t 25]. 

6 

The angle 0 between the B meson momentum direction and the beam axis in the CM frame. 

7 

The angle between the D and T(4A) directions in the rest frame of the B candidate. 

8 

The projection of the sphericity vector with the largest eigenvalue onto the e + e - beam direction. 

9 

The angle of the sphericity vector with the largest eigenvalue with respect to that of the remaining 

particles. 

10 

The angle of the sphericity vector with the second largest eigenvalue. 

11 

The angle of the sphericity vector with the smallest eigenvalue. 

12 

The magnitude of the thrust of the particles not used to reconstruct the signal. 


TABLE I: Variables used for qq suppression. 


removes 19% of cc background and 11% of BB background according to MC simulation. The relative loss in signal 
efficiency is 5.5%. 

In the rare case where there are multiple candidates in an event, the candidate with Mb c closest to the nominal 
value is chosen. The relative loss in signal efficiency is 0.8%. 

To discriminate signal events from the large combinatorial background dominated by the two-jet-like e + e~ -A qq 
continuum process, where q indicates u, d , s or c, a multivariate analysis is performed using the 12 variables introduced 
in Table [TJ To effectively combine these 12 variables, we employ the NeuroBayes neural network package [26] . The 
NeuroBayes output is denoted as Cnb and lies within the range [—1, 1]; events with Cnb ~ 1 are signal-like and 
events with Cnb ~ — 1 are (/(/-like. Training of the neural network is performed using signal and qq MC samples. 
The Cnb distribution of signal events peaks at Cnb ^ 1 and is therefore difficult to represent with a simple analytic 
function. However, the transformed variable 


Cnb — l n 


Cnb ~ Cnb ,low 
Cnb ,high Cnb 


(9) 


where Cnb, low = —0.6 and Cnb, high = 0.9992, has a distribution that can be modeled by a Gaussian for signal as well 
as background. The events with Cnb < Cnb, low are rejected; the relative loss in signal efficiency is 7.4%. 


ANALYSIS PROCEDURE 

In this section we describe the fit to determine the physics parameters. In Sec. [] we describe the signal and 
background shape parametrization. In Sec. [] we describe how we correct for the effect of migration and acceptance 
variations between bins. In Sec. [] the fit to extract the values of (x, y) is described. 


Signal and background parametrization 

The number of signal events is obtained by fitting the three-dimensional distribution of variables Mb c , A E, and 
C^ B using the extended maximum likelihood method. We form three-dimensional PDFs for each component as the 
product of one-dimensional PDFs for A E, M^ c and C^ B , since the correlations among the variables are found to be 
small. The fit region is defined as A E G [—0.1,0.15] GeV and Mb c > 5.21 GeV/c 2 . 

Backgrounds are divided into the following components: 

• Continuum background from qq events. 

• BB background, in which the tracks forming the B° -a DK*° candidate come from decays of both B mesons in 
the event. The number of possible B decay combinations that contribute to this background is large; therefore, 
both the Dalitz distribution and distribution of the fit parameters are quite smooth. BB backgrounds are 
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FIG. 3: Projection of the fit to real data using the full Dalitz plot. Left: A E distribution with Mb c > 5.27 GeV/c 2 and 
C^b > 2 requirements. Middle: C^b distribution with | A E |< 0.03 GeV and Mb c > 5.27 GeV/c 2 requirements. Right: Mb c 
distribution with | A E |< 0.03 GeV and C^b > 2 requirements. Curves show the fitted signal and background contributions, 
(red is signal, yellow is D°af , green is D°p° , blue is Df a ke BB, light blue is D true BB and magenta is qq) and points with error 
bars are the data. 


further subdivided into two components: events reconstructed with a true D decay, referred to 

as Dtrue BB background, and those reconstructed with a combinatorial D candidate, referred to as Df a k e BB 
background. 

• Peaking BB background, in which all tracks forming the B° —»• DiG*° candidate arise from the same B meson. 
This background has two types: events with one pion misidentified as a kaon, such as D°[ 7 r + 7 r _ ] p o, and one pion 
misidentified as a kaon and one pion not reconstructed, such as D°[ 7 r + 7 r + 7 r _ ] a +. The backgrounds come from 
individual B decays and are well separated from the signal. 

The A E PDFs are parameterized by a double Gaussian for the signal, an exponential function for the D true BB 
background, an exponential function for the Df a k e BB background, a linear function for the qq background, a double 
Gaussian for the D°p° background, and a Gaussian for the D°af background. The M^ c PDFs are a Gaussian for 
signal, a Crystal Ball function [27. for D true BB background, an ARGUS function [28| for Df a k e BB background, an 
ARGUS function for qq background, a sum of a Gaussian and ARGUS functions for D° p° background and a Gaussian 
for D°al background. For each component, the C^ B PDF is the sum of a Gaussian and bifurcated Gaussian. The 
shape parameters of the PDFs are fixed from MC samples. 

The numbers of events in each bin are free parameters in the fit. This procedure has been justified for background 
that is either well separated from the signal (such as peaking BB background) or is constrained by a much larger 
number of events than the signal (such as qq background). The results of the fit to the full Dalitz plot are shown 
in Fig. [ 3 J We obtain a total of 44.2^^ 1 signal events. The statistical significance is 2 .8a relative to the no-signal 
hypothesis. Simultaneously, we obtain 695.8l^5 g for D true BB, 1963.2^227!5 f° r ^fake BB, 11075. 7 ^ 55;5 for qq , 
16.6^3*6 for D°p° and 59.3^20 8 f° r backgrounds events. 


Corrections to the bin-by-bin yields 

There are further effects that must be accounted for before the values of (x, y ) can be determined from a binned fit. 
Equation [7] only holds if there is no migration between bins and the Dalitz acceptance is uniform. Here, we consider 
the crossfeed for bin-by-bin yields as the migration and the acceptance as the event reconstruction efficiency. 

First, we discuss migration which is due to momentum resolution and flavor misidentification. Momentum resolution 
leads to migration of events among the bins. In the binned approach, this effect can be corrected in a non-parametric 
way. The migration can be described by a linear transformation of the number of events in each bin 

N 0 bs,i = y^ikNj ( 10 ) 

where N- is the the number of events that bin i would contain without the migration with acceptance and 7V 0 b s ,i is 
the reconstructed number of events in bin i. The migration matrix is nearly the unit matrix; it is obtained from 






















































a signal MC simulation generated with the amplitude model reported in Ref. [5]. Most of the off-diagonal elements 
are null; only a few have values \a^\ < 0.04. In the case of a D decay from a £>, the migration depends 

on the parameters x and y. However, this is a minor correction to an already small effect and so is neglected. 

The second migration effect to be considered is due to misidentification of the B flavor. Double misidentification 
in K*° reconstruction from K + tt~ where K~ is misidentified as tt~ and 7r + is misidentified as K + at the same time 
leads to migration of events between N+ NZi due to assignment of the wrong flavor to the B candidate. If the 
fraction of doubly-misidentified events is /?, the number of events in each bin can be written as 

^ = + ( n ) 

The value of /3 is obtained from MC simulation and is found to be (0.12 d= 0.01)%. Therefore, the effect of flavor 
misidentification is neglected. 

The final state radiation also causes migration between bins. The measured values of c* and Si by CLEO are not 
corrected for the radiation and the effect upon our analysis is found to be negligible m ■ 

Second, we consider the effect of the variation of the efficiency profile over the Dalitz plane. We note that Eq. © 
does not change under the transformation P —>> eP when the efficiency profile e(m\,rn 2 _) is symmetric: = 

The effect of non-uniform efficiency over the Dalitz plane cancels when using a flavor-tagged D sample 
with kinematic properties that are similar to the sample from the signal B decay. This approach allows for the removal 
of the systematic uncertainty associated with the possible inaccuracy of the detector acceptance description in the 
MC simulation. With the efficiency taken into account (that is, in general non-uniform across the bin region), the 
number of events reconstructed is 

N' = j p{V)e(V) dV. (12) 

Here, p is the probability density on the Dalitz plane and V is the position on Dalitz plane. Clearly, the efficiency 
does not factorize. One can use an efficiency averaged over the bin, then correct for it in the analysis: 


N[ _ f p(V)e(V)dV 

N~i ~ JpCDjdD 


(13) 


Here, W are the number of events corrected for variations in acceptance and migration, which should be used for (x±, 
y±) extraction. The averaged efficiency e* can be determined from MC. The assumption that the efficiency profile 
depends only on the D momentum is tested using MC simulation and the residual difference is treated as a systematic 
uncertainty. The correction for c$ and Si due to efficiency variation within a bin cannot be calculated in a completely 
model-independent way, since the correction terms include the amplitude variation inside the bin. Calculations using 
the Belle D —>> model [5] show that this correction is negligible even for very large non-uniformity of the 

efficiency profile. 


Fit to determine (x , y) 

If Ni in each bin is measured, x± and y± can be obtained according to Eq. 0 by minimizing 

- 21og£(x,y) = -2^1ogp((ATi) ( x,y),N h a Ni ), (14) 

where (Ni) are the expected number of signal events in the bin i obtained from Eq. 0 . Here, Ni and a n. are the 
observed number of events in data and the uncertainty on W, respectively. 

The procedure described above does not make any assumptions about the Dalitz distribution of the background 
events, since the fits in each bin are independent. Thus, there is no uncertainty related to the Dalitz model. However, 
in our case, where there are a small number of events and many background components, such independent fits are 
not feasible. Therefore, we obtain (x±, y±) from a combined fit with a common likelihood for all bins. The relative 
numbers of background events in each bin are constrained to the numbers found in the MC. The amount of the 
D true BB background in bins from the ratio of D° ( Ki ) and D° (K-i) from MC and the amount of the Df a k e , BB , 
qq and the background from individual B decays from the MC. The yields integrated over the Dalitz plot of the 
background components are additional free parameters. Thus, the variables (x±, y±) become free parameters of the 
combined likelihood fit and the assumption that the signal yield obeys a Gaussian distribution is not needed. While 
the normalization parameter Kb is also a free parameter of the fit, we do not mention it in the following as it is not 
a quantity of interest. 
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FIG. 4: Projections of the combined fit of the B° —> DK*° sample on A E for each Dalitz bin, with the M}- IC > 5.27 GeV/c 2 
and Cnb > 2 requirements. The fill styles for the signal and background components are the same as in Fig. [ 3 ] 


COMBINED FITS TO DATA 

The results of the combined fit in each bin of the B° and B° are shown in Figs. [ 4 ] and [ 5 J respectively. The 
plots show the projections of the data and the fitting model on the A E variable, with the additional requirements 
Mbc > 5.27 GeV/c 2 and C^ B > 2. The values of the (x, y ) parameters and their statistical correlations, obtained 
from the combined fit for the signal sample, are given in Table |II| In this study, these ( x , y) values from the likelihood 
distribution of the combined fit are corrected using the frequentist approach with Feldman-Cousins ordering [29 ; , 
which is described in Sec. fl 
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sample on AE for each Dalitz bin, with the Mb c > 5.27 GeV/c 2 
and > 2 requirements. The fill styles for the signal and background components are the same as in Fig. [ 3 ] 


SYSTEMATIC UNCERTAINTIES 


The systematic uncertainties of (x, y) are obtained by taking deviation from the default procedure under various 


assumptions. The systematic uncertainties are summarized in Table III most are negligible compared to the statistical 
uncertainty. There is an uncertainty due to the Dalitz efficiency variation because of the difference in average efficiency 
over each bin for the flavor-tagged D and B° —>■ DK*° samples. A maximum difference of 1.5% is obtained in a MC 
study. The uncertainty is taken as the maximum of two quantities: 


• the root mean square of x and y from smearing the numbers of events in the flavor-tagged sample K{ by 1.5%, 
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Parameter 


X- 

+0.29 ±0.32 

y- 

-0.33 ±0.41 

corr.(x, y-) 

±7.0% 

x+ 

±0.07 ±0.42 

y+ 

±0.05 ± 0.45 

corr.(x+, y + ) 

-7.5% 


TABLE II: (x, y ) parameters and their statistical correlations from the combined fit of the B° —» DK*° sample. The error is 
statistical. The values and errors are obtained from the likelihood distribution. 


Source of uncertainty 

Ax- 

Ay- 

Ax+ 

^2/+ 

Dalitz efficiency 

±0.00 

+ 0.01 

- 0.00 

±0.01 

+ 0.00 

- 0.01 

Migration between bins 

±0.00 

+ 0.01 

- 0.00 

+ 0.01 

- 0.00 

±0.00 

PDF parameterization 

+ 0.01 

- 0.07 

+ 0.07 

- 0.01 

+ 0.01 

- 0.10 

+ 0.04 

- 0.06 

Flavor-tag statistics 

±0.00 

±0.00 

±0.00 

+ 0.00 

- 0.01 

Ci, Si precision 

±0.03 

+ 0.09 

- 0.08 

±0.05 

+ 0.08 

- 0.10 

k precision 

±0.00 

±0.01 

±0.00 

±0.00 

Total without a, Si precision 

+ 0.01 

- 0.07 

+ 0.07 

- 0.02 

+ 0.02 

- 0.10 

+ 0.04 

- 0.06 

Total 

+ 0.03 

- 0.08 

+ 0.12 

- 0.08 

+ 0.05 

- 0.11 

+ 0.09 

- 0.12 


TABLE III: Systematic uncertainties in the (x, y ) measurement for the B° -+ DK*° mode. Values are rounded to two significant 
digits and those less than 0.005 are quoted as 0.00. 


or 

• the bias in x and y between the fits with and without efficiency correction for A+ obtained from signal MC. 

The uncertainty due to migration of events between bins is estimated by taking the bias between the fits with and 
without the migration correction. The uncertainties due to the fixed parameterization of the signal and background 
PDFs are estimated by varying them by ±1 a. The uncertainty due to the C^ B PDF distributions for BB is estimated 
by replacing them with the signal C^ B PDF. The uncertainty due to the D true and D^ e BB fractions is estimated 
by varying them between 0 and 1. The uncertainty arising from the finite sample of flavor-tagged D -+ K$tttt decays 
is evaluated by varying the values of A+ within their statistical uncertainties. The uncertainty due to the limited 
precision of <+ and s* parameters is obtained by smearing the c* and Si values within their total errors and repeating 
the fits for the same experimental data. The uncertainty due to k in Eq. © is evaluated by varying the value of 
k(= 0.95 ± 0.03) within its error [16]. Total systematic uncertainties in the (x, y) are obtained by summing all 
uncertainties in quadrature and listed in Table m 


RESULT 


We use the frequentist approach with Feldman-Cousins ordering [29] to obtain the physical parameters y = 
( 03 ,rs , Ss) (or true parameters fi = z trU e = y~, £+,?/+)) from the measured parameters z = z meas = 

(x_, x + j ?/_(_) taken from the likelihood distribution. In essence, the confidence level a for a set of physical param¬ 
eters y is calculated as 

, s _ fv( M )P( Z I P) dz 
aM I 00 p{z\fj)dz 

where p(z \ y) is the probability density to obtain the measurement z given by the set of true parameters y. The 
integration domain V(y) is given by the likelihood ratio (Feldman-Cousins) ordering: 


(15) 


P(z I M) 

P(z | Aibest(^)) 


> 


pjzp I jj) 

P(zo I Mbest(^o))’ 


( 16 ) 
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FIG. 6: C.L. contours for (x_, y~) (blue) and (x+, y+) (red). The dots show the most probable (x, y) values; the lines show 
the 68% contours. The fluctuations arise from the statistics of the pseudo-experiments and C.L. step used. 


where p best ( 2 ) is p that maximizes p(z \ p) for the given z, and zq is the result of the data fit. This PDF is taken 
from MC pseudo-experiments. 

Systematic uncertainties in p are obtained by varying the measured parameters z within their systematic uncer¬ 
tainties assuming a nominal distribution. In this calculation, we ignore the correlations of uncertainties between the 
B° and B° as the two samples are independent. 

As a result of this procedure, we obtain the confidence levels (C.L.) for (x, y) and the physical parameter rs- The 
C.L. contours on (x, y) are shown in Fig. [6j The 1 — C.L. as a function of rs is shown in Fig. [ 7 ] The final results are: 


X- 

y- 

x+ 

y+ 

r s 


< 


+0.4 


+1.0+0.0 
-0.6-0.1 


± 0 . 0 , 


- 0.6 


+0.8+0.1 
-1.0-0.0 


± 0 . 1 , 


+0.ltg;Itg;S±0.l, 

+0.3tg;iig;S±0.1, 

0.87 at 68% C.L., 


(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


where the first error is statistical, the second is systematic without uncertainties in (c$, s^), and the third is from the 
(ci, Si) precision from CLEO. 


CONCLUSION 

We report the first measurement of the amplitude ratio rs using a model-independent Dalitz analysis of D -+ 
Ks7t^tt~ decays in the process B° -+ DK*° with the full data sample of 711 fb -1 corresponding to 772 x 10 6 BB 
pairs collected by the Belle detector at the T(4*S) resonance. Model independence is achieved by binning the Dalitz 
plot of the D —>> iFg7r + 7r - decay and using the strong-phase coefficients with binning as in the CLEO experiment m- 
We obtain the value rs < 0.87 at 68% C.L. This measurement results in lower statistical precision than the model- 
dependent measurement from BaBar with the B° -+ DK° mode [8] despite the larger data sample due to the smaller 
B° -+ DiF*° signal observed. The result is consistent with the most precise rs measurement reported by the LHCb 
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FIG. 7: Likelihood profile for rs- The blue points are for B° (x_, y_), red are for B° ( 2 ;+, y+) and black are B° and B° 
combined. The two horizontal lines show 68% and 95% C.L. 


Collaboration [17] of rs = 0.2401q!o 48 that uses [K + K~, K ± k t , 7 t^7t~]dK*° decays. We have confirmed 

the feasibility of the model-independent Dalitz analysis method with neutral B —)► DK*. The value of rs indicates 
the sensitivity of the neutral B —» DK* decay to (j )3 because the statistical uncertainty is proportional to 1 /rs- In 
future high statistics experiments such as Belle II and the LHCb upgrade, this method will give a precise and model 
independent determination of </> 3 . A more advanced double-Dalitz-plot analysis of B° DK + 7 r _ , D —>> Ks7r^7r~ - ED] 
has been proposed; this result can be considered as that from one bin of such an analysis. 
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